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ABSTRACT 

Using CHARA and VLTI ncar-infrarcd spcctro-intcrfcromctry with hcctomctric baseline lengths (up 
to 330m) and with high spectral resolution (up to A/AA = 12 000), we studied the gas distribution 
and kinematics around two classical Be stars. The combination of high spatial and spectral resolution 
achieved allows us to constrain the gas velocity field on scales of a few stellar radii and to obtain, for 
the first time in optical intcrferometry, a dynamical mass estimate using the position- velocity analysis 
technique known from radio astronomy. For our first target star, (3 Canis Minoris, we model the 
H-f K-band continuum and Br7-linc geometry with a near-critical rotating stellar photosphere and 
a geometrically thin equatorial disk. Testing different disk rotation laws, we find that the disk is 
in Keplerian rotation {v{r) cx 7.-0-5±o.ij ^^^^ derive the disk position angle (140 ± 1.7°) inclination 
(38.5 ± 1°), and the mass of the central star (3.5 ± 0.2 Mq). As a second target star, we observed the 
prototypical Be star C Tauri and spatially resolved the Br7 emission as well as nine transitions from 
the hydrogen Pfund series (Pf 14-22). Comparing the spatial origin of the different line transitions, 
wc find that the Brackett (Br7), Pfund (Pfl4-17), and Balmcr (Ha) lines originate from different 
stellocentric radii (i?cont < ^Pf < Rbi-j Riia), which we can reproduce with an LTE line radiative 
transfer computation. Discussing different disk-formation scenarios, we conclude that our constraints 
are inconsistent with wind compression models predicting a strong outflowing velocity component, but 
support viscous decretion disk models, where the Keplerian-rotating disk is replenished with material 
from the near-critical rotating star. 

Subject headings: circumstellar matter - stars: emission-line, Be - stars: individual (/3 CMi. ^Tau) — 
stars: fundamental parameters - techniques: interferometric 



1. INTRODUCTION 

Classical Be stars are main-sequence (or near main- 
sequence) B-type stars associated with hydrogen line 
emission, indicating the presence of ionized circumstel- 
lar gas, which is believed to be arranged in an equato- 
rial disk- like structure. Optical/infrared spectroscopic 
observations have r evealed non-radial pulsations (e.g. 
IRivinius et al.|[l998[ ). which might provide a way to feed 
material from the ph otosphere to the in ner disk. Po- 
larimctric studies (e.g. iDraper et al.ll20lTI ) can constrain 
the disk density structure. Further unique insights into 
the structure and physics of these disks can also be ob- 
tained with intcrferometry at visual and infrared wave- 
lengths, allowing one to unravel the inner disk structure 
on scales of a few stellar radii directly. For instance, 
interferometric studies have allowed associating quasi- 
cyclic variations in the ratio between the blue- and red- 
shifted wing of the Ha-linc emission (V/R variability) 
with global oscillations in the circumstellar disk, likely 
in the form of a one-armed spiral density pattern (e.g. 
IVakiH et al.l[l99l: IStefl et al.l[2009t ICarciofi et al.l[200l 



iSchaefer et al.|[2010[ ). Besides studies on the disk con- 
tinuum geometry, interferometric observations in spec- 
tral lines have provided the first direct constraints on 
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^ Based on observations made with ESO telescopes at the 
Paranal Observatory under programme IDs 084.C-0848(A) and 
085.C-0911(A) and with the CHARA array 



evidence that the disks around classical Be stars exhibit 
a near-Keplerian rotation profile, which might allow to 
effectively ru le out several disk-formation scenarios (e.g. 
see review bv lCarciofill2010( ). However, most earlier stud- 
ies using spectro-interferometry were limited in terms of 
spectral resolution or baseline position angle (PA) cover- 
age, leaving significant uncertainties about the detailed 
gas velocity field and the evidence to distinguish between 
a purely rotational versus an expanding velocity com- 
ponent in the disk. Obtaining such evidence is essen- 
tial in order to decide between different scenarios which 
have been proposed to explain the disk formation mecha- 
nism, including radiatively driven winds, ram pressure or 
magneti cally induced wind comp ression, and viscous de- 
cretion (jPorter fc Riviniuil2003f l. Furthermore, recently 
there has been a controversy about the appearance of 
a phase inversion in spectro-interferometric observations 
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of several classical Be stars, which triggered speculations 
about secondary dynamical effects or the need for an 
additional kinematic al component bey ond the canonical 
star+disk paradigm (jStefi et alJl201l( i. 

Here, wc present near-infrared spectro-intcrfcromctric 
observations on the classical Be star /3 CMi with a high 
spectral resolution oi R = 12, 000 in the hydrogen Br7- 
line, enabling us to constrain the rotation profile directly. 
In addition, we observed the classical Be star (^Tau for 
the first time in multiple hydrogen line transitions (Br7 
and Pfund lines), providing direct information about the 
excitation structure within the disk. 

In the following, we present our CHARA and VLTI 
interferometric observations (Sect. [2]). The observations 
in spectral lines are then first interpreted using a model- 
independent photocenter analysis approach (Scct.|3]). In 
Sect.m we present our continuum and kinematical mod- 
eling on /3 CMi, followed by our discussion of the results 
on (^Tau (Sect. [5]). Finally, we summarize our findings 
in Sect. [6l 

2. OBSERVATIONS 
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Fig. 1. — Jiti-coverage achieved with our CHARA/MIRC {H- 
band), CHARA/CLIMB (it-band), and VLTI/AMBER (it-band) 
interferometric observations on /3 CMi. 

Our near-infrared _ff-band continuum observations 
on /3 CMi were obtained u sing the CHARA array 
(|ten Brummelaar et al.l I2005D . which is operated by 
Georgia State Univers ity. The MIRC beam combiner 
onnier et al.ll2006b[ ) allowed us to combine the light 
from four of the six CHARA 1 m-telescopcs simultane- 
ously, yielding baseline lengths of up to 330 m, with good 
baseline coverage (Fig. [T]). The MIRC data cover the 
77-band with low spectral dispersion {R = 35) and was 
reduced using th e University of Michig an MIRC data re- 
duction pipeline (jMonnier et al.l[2007[ ). 

In order to investigate the /iT-band disk geometry of 
/3CMi, we empl oyed the CHARA/CL IMB 3-telescope 
beam combiner (jSturmann et al.ll2010D . Visibilities and 
closure phases were derived using the "redclimb" and 
"reduceir" software. Besides the statistical errors, we 
also add a calibration uncertainty of 0.05 for the derived 
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Gas distribution, kinematics, and excitation structure in the disks around /? CMi and C Tau 
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Fig. 2. — Closure phases measured with the MIRC beam com- 
biner on /3 CMi (blue data points), overplottcd with the model 
predi ction s from our best-fit _ff-band model (red data points; 
Sect. 11211. 

visibilities, which represents an empirical value for the 
typical scatter in the instrument transfer function. 

Spectro-interferometric observations with medium 
(MR mode, R = 1500) and high spectral dispersion (HR 
mode, R ~ 12 000) were obtained with the Very Large 
Telescope Interferometer (VLTI) of the European South- 
ern Observatory a nd the AMBER 3-t elescope beam com- 
biner instrument (jPetrov et al.ll2007[ ). The AMBER ob- 
servations on /3 CMi were recorded using three 8.2 m 
unit telescopes (2009-12-31) and three 1.8 m auxiliary 
telescopes (2010-04-23), respectively. The atmospheric 
piston was stabilized us ing the FINITO fringe tracker 
(|Le Bouquin et al1l2008| ). which allowed us to use long 
detector integration times (DITs) of 1 s and 6 s and to 
record data with high spectral dispersion around the 
hydrogen Br7-line (A™™""^ = 2.166078 urn). Due to 
the presence of residual phase jitter, the absolute vis- 
ibility calibration of the AMBER data is not reliable, 
while the important wavelength-differential observables 
are not affected. Spectra and wavelength-differential vis- 
ibilities and phases (DPs, Fig. [3]) were derived from the 
AMBER data using the am d lib (V3.0) data re duction 
software (|Tatulli et all [20071: iCheUi et all [2009I ). The 
wavelength calibration was done using atmospheric tel- 
luric features close to the Br7-line (yielding an accu- 
racy of approximately 1 spectral channel) and by ap- 
plying a heliocentric-barycentric system correction us- 
ing heliocentric velocities of 4-5.97 kms""'^ (2009-12-31) 
and -28.61 kms'^ (2010-04-23), respectiv ely For the 
syste mic velocity we assume +22.0 kms-i (jPuflot et al.l 
[19951) . 

From all CHARA and VLTI interferometric observa- 
tions, we also derived closure phases (CPs). Both the 
H- and K-hand continuum closure phases are consistent 
with zero on a 2(T-level, which leads us to conclude that 
the brightness distribution does not show significant indi- 
cations for deviations from centro-symmetry. The most 
constraining CPs have been recorded with MIRC in the 
iJ-band, which are shown in Fig. [2 Our AMBER HR 
measurements from 2009-12-31 provide us also with a CP 
measurement in the Br7-line of /? CMi, while the derived 
CPs from the 2010-04-23 dataset are rather noisy and 
are therefore not included for our model fits. 

CTau was observed on 2010-01-01 using AMBER's 



MR-mode covering the upper A'-band (2.12 to 2.46 fim). 
The assumed syste mic velocity for CTau is +21.8 kms^^ 
(jPuflot et al.lll99"5l ). In the ^Tau data, we detected not 
only the Br7 transition (7-4), but also hydrogen Pfund 
transitions (Fig. |4l 1st row), including clear detection 
of Pfl4 (19-5, 2.4477 /itm) to Pf22 (27-5, 2.3591 fim). 
Higher Pfund transitions are also present in the spec- 
trum, but cannot be clearly separated due to the wide, 
double-peaked profile of the individual lines. For C Tau, 
the derived CPs are too noisy to provide additional in- 
formation and are therefore not included in our further 
analysis. 

Details about the observational setup for all interfero- 
metric observations are listed in Tab.[T] Each observation 
on a science star was accompanied by observations on in- 
terferometric calibrators, allowing us to monitor and cor- 
rect for the atmospheric and instrumental transfer func- 
tion. 

3. PHOTOCENTER ANALYSIS 

Differential phases provide unique information about 
small-scale (sub-mas) photocenter displacements be- 
tween the blue- and red-shifted line wings. These dis- 
placements provide a very sensitive measure of the gas 
kinematics on scales of a few stellar radii. In a first anal- 
ysis step, we reconstruct the on-sky 2-D photocenter dis- 
placement from the measured DPs by solving the system 
of linear equations 

^ (t)i X 

where is the differential phase measured on baseline 
i, Bi is the correspo nding baseline vector, a nd A is the 
central wavelength (|Le BouQuin et al] |2009() . The de- 
rived photocenter plots for /3 CMi and C Tau are shown in 
Figs.[3]and|4|(mid(i/e panel), respectively, and clearly re- 
veal rotation-dominated velocity fields for both objects, 
as indicated by the linear alignment of the photocenter 
vectors for the different gas velocities and the opposite 
sign of the photocenter displacement for the blue- and 
red-shifted emission. In order to associate the DP with 
the on-sky orientation, it is necessary to calibrate the 
sign of the DP measurements. For this purpose, we re- 
processed the CTau VLTI/ AMBER data set presented 
bv lStefl et all (|2009l) and calibrated our DP sign in order 
to match the published on-sky orientation. 

For (3 CMi, we determine the position angles 9 for the 
disk rotation plane at the two epochs to be 138.3 ± 1.5° 
(2009-12-31) and 141.7+1.4° (2010-04-22). The photo- 
center vectors corresponding to our highest-SNR obser- 
vation (2009-12-31) show an interesting arc-like structure 
in the red-shifted line wing (Fig. [3l bottom left), where 
the photocenter vectors corresponding to low gas veloci- 
ties are above the derived disk plane, while the high ve- 
locities are displaced in the opposite direction. Although 
the significance of this pattern is still only marginal in 
our data, we speculate that this pattern might result 
from opacity effects, with the more distant parts of the 
disk appearing fainter than the disk parts facing the ob- 
server. Such an obscuration screen would displace the 
photocenter perpendicular to the disk plane, where the 
amplitude of the displacement is stronger for lower gas 
velocities, since the low-velocity emission is distributed 
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Fig. 3. — Upper panel: VLTI/AMBER spectra (1st row), visibilities (2nd row), and DPs (3rd row) measured on /3 CMi for tlie epochs 
2009-12-31 (left) and 2010-04-22 (right). Middle panel: From the measured DPs, we derive for each spectral channel the 2-D photocenter 
displacement vector (including continuum and line contributions; East is plotted left and North is up). The flux and DP data points as 
well as the derived photocenter vectors have been color-coded based on the Doppler velocity (see the upper panel to relate each color to a 
wavelength). The DPs corresponding to the astrometric solutions are shown in the middle panel as solid black lines. Bottom panel: Using 
the procedure outlined in Sect. \3\ wc corrected for the continuum contributions in the line spectral channels, revealing the photocenter 
displacement corresponding to the line emission only. Besides the determin ed p hotocenter offsets, for comparison we also show th e siz e of 
the //-band continuum-emitting disk (black ellipse, FWHM Gaussian, Sect. l4?2ll and of the stellar photosphere (grey ellipse. Sect. |4!T1| . 
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Fig. 4. — Upper panel: VLTI/AMBER spectra (1st row), visibilities (2nd row), and DPs (3rd row) measured on (^Tau for the epoch 
2010-01-01 in the Br7 (left) and Pfl4-22 transitions (right). Middle/bottom panel: Photocenter displacement vectors derived from the 
measured DPs (including line and continuum emission, middle panel) and the continuum-corrected DPs (traci ng the line emi s sion o nly, 
bottom panel). For comparison, we also show the size of the //-band continuum-emitting disk, as determined bv lSchaefer et all I I2010I) for 
epoch 2009-11-10 (black ellipse, FWHM Gaussian) and of the stellar photosphere (grey ellipse). 
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over a more extended region. Accordingly, the displace- 
ment would be strongest at zero velocities, and then sym- 
metrically decrease towards higher velocities. The super- 
position of this weak displacement (perpendicular to the 
disk plane) with the displacement due to Kcplcrian ro- 
tation (parallel to the disk plane) might result in the 
observed arc-shaped structure, which would also provide 
a unique tool to determine the orientation of the disk in 
space and the disk rotation sense. In the case of /3 CMi, 
this implies that the north-eastern part of the disk is fac- 
ing towards the observer (based on the displacement of 
the low-velocity channels in this direction) and that the 
disk is in clockwise rotation (based on the location of 
the red-shifted photoccnter displacements in the north- 
western quadrant). 

Another intriguing feature in the /3 CMi data is the 
M^-shaped profile, which we observe in the wavelength- 
dependent visibilities and DP on our longest interfero- 
metric baselines (Fig. [31 2nd and 3rd row). Likely, this 
profile indicates that the visibility function of the line- 
emitting region passes through a visibility null and tran- 
sits from the first to the second visibility lobe. Since 
this effect would reverse the direction of the photoccn- 
ter vector, we manually correct the phase sign in these 
corresponding spectral channels close to the line center. 
For C Tau. we can determine the rotation axis for the 
Br7 (123.9 ± 2.0°) and the nine Pfl4-Pf22 transitions 
(127.6 ± 1.5°) separately and find that the line-emitting 
gas rotates in the same disk plane within the observa- 
tional uncertainties of ^ 2°. 

Using the aforementioned procedure, it is possible to 
reliably measure the direction of the photocenter dis- 
placement, while the length of the displacement vector 
is biased by the contributions from the underlying con- 
tinuum emission. In order to remove these contributions 
from the measured observables (F, V, ( /)) , we apply the 
method outlined bv lWeigelt et al.l (|2007[ ) and interpolate 
the continuum flux (Fc) and continuum visibility (Fc) 
from the adjacent continuum. The visibility and DP of 
the pure line emitting- region (VJ, (/>;) are then given by 

\FiVi\''^\FV\^ + \F,Vc\^ ~2-FV- F,V^ ■ cos<^ (2) 
\FV\ 

sin0i=sin(?!)|^^, (3) 

where Fi = F — Fc denotes the flux contribution from 
the spectral line. The continuum-corrected DPs arc then 
used to derive the photocenter displacement of the pure 
line-emitting region (Figs. 13] and 01 bottom panel). Ap- 
plying this correction will provide the real centroid off- 
set of the line emission, but also introduce noise from 
the visibility and flux measurements, resulting in an in- 
creased scatter in the position angle distribution. There- 
fore, we decided to measure the position angle of the ro- 
tation axis from the uncorrected line-l-continuum photo- 
center displacements (Figs. |31 and SI middle panel), while 
the continuum-corrected photocenter displacements (bot- 
tom panel) will be used later on to construct a position- 
velocity diagram and to compare the spatial origin in 
different line tracers. 

4. DISCUSSION ON 13 CANIS MINORIS 
/3CMi (HR2845) is a relatively quiet B8V-type clas- 
sical Be star located at a distance of 52.2I2.2 PC 



(|Tvcner et al.l [20051 ). Ground-based photometric moni- 
toring has provided no clear ind ications for signiflcant 
variability (jPavlovski et al.l I1997D , while there is some 
margin al evidence for long-term variati ons in the Ha- 
profilc ('PoUm annll2002l : lHesselbachll2009( ) . 

In the following, we present a refined model for the 
photospheric emission of (3 CMi taking the near-critical 
stellar rotation into account (Sect. ITT]) , followed by our 
modeling of the interferometric data in continuum emis- 
sion (Sect. |4?2|) . Given that the disk rotation signatures 
in the Br7 emission line are overlayed on the rotation 
signatures of the star (in Br7 absorption) , we investigate 
the influence of the stellar rotation on our measurements 
(Sect. 14. 3p . In the following, we construct a position- 
velocity diagram (Sect. 14. 4|) and present a full kinemati- 
cal modeling in Sect. 14.51 

4.1. Constraining the stellar parameters 

In order to obtain a model for the photospheric emis- 
sion of /3 CMi, we searched for a consistent set of stellar 
parameters, taki ng the evidence for near-critical rota- 
tion into account (jSaio et al.ll2007D . For this purpose, we 
employed our r apid rotator code (jMonnier et al.l [20071 : 
iChe et aI1l2011h . which simulates the stellar oblateness 
and surface temperature distribution using the modi- 
fied von Zeipel theorem (gravity darkening coefficient = 
0.188), and computes model images in the continuum 
and in the photospheric Br7 absorption line. 

In order to constrain the stellar parameters, we com- 
puted a small parameter grid, in which we systemati- 
cally varied the inclination angle i and the fractional 
angular velocity w/wcrit (where Wcrit is the critical an- 
gular velocity). For a given inclination angle and frac- 
tional angular velocity, we start our iterative process 
by assuming a stellar mass M^, from which we com- 
pute the polar radius i?poio using the rotation velocity 
vsini = 244 ±6 kms'^ (|Yudinll200i[ l. The polar tem- 
perature Tpoic is derived from the y-band magnitude 
{V ~ 2.89). Any disk continuum flux in the optical 
band will tend to dilute the photospheric spectrum caus- 
ing absorption lines to appear weaker than expected. 
We checked for such line reduction in blue spectra of 




4460 4470 4480 
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Fig. 5. — Spectrum of the He I 447.1 nm and Mg II 448.1 nm line 
of /3 CMi (solid line) . The dashed line shows a photospheric model 
without disk contributions, while the dotted line corresponds to a 
model where 20% disk continuum contributions have been added. 
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TABLE 2 

Grid of stellar parameters, obtained for /3 CM: using our fast rotator model (Sect. [iTTjl . 



Model 


i 

[°] 




[A/q] 


-^polc 

[Re] 


[Rq] 


^polc 

[K] 


P 

[days] 


[Le] 


[K] 


A 


45 


0.99 


3.70 


3.32 


4.61 


13200 


0.676 


267 


11558 


B 


45 


0.96 


4.11 


3.03 


3.94 


14600 


0.577 


343 


13085 


C 


45 


0.93 


4.47 


2.87 


3.59 


15650 


0.527 


413 


14235 


D 


40 


0.99 


3.98 


2.95 


4.10 


14550 


0.546 


337 


12871 


E 


40 


0.96 


4.49 


2.74 


3.56 


16050 


0.475 


435 


14505 


F 


40 


0.93 


4.91 


2.60 


3.25 


17300 


0.433 


532 


15847 


G 


35 


0.99 


4.43 


2.62 


3.64 


16200 


0.433 


437 


14447 


H 


35 


0.96 


5.03 


2.44 


3.17 


18100 


0.377 


590 


16474 



pCMi obtained by iGrundstroml (|2007| ). Fig. [5] shows 
the observed profiles {solid line) of He I 447.1 nm and 
JVIg II 448.1 nm based upon the average of 11 spectra 
made between 2005 and 2009 with the Kitt Peak Na- 
tional Observatory Coude Feed Telescope (spectral reso- 
lution R = 12, 500). The observed spectrum is compared 
with a theoretical flux spectrum (dashed line) from th e 
grid of LTE models by lRodrfguez-]V[erino et al.l ()2005D . 
We adopted average parameters for the visible hemi- 
sphere of the star of Tpfp = 11,800 K, l ogq = 3.8, and 
y sin z = 230 km s~ 1 (jFr^mat et al.ll2005n , and the model 
spectrum was convolved with a simple rotational broad- 
ening functio n for a line ar limb darkening coefficient of 
e = 0.42 (jWade fc Rucinski 19851 ). Note that this ap- 
proach assumes a spherical star, ignores gravity darken- 
ing, and neglects changes in the local intensity spectrum 
with orientation, but these simplifications are acceptable 
for our purpose of checking for systematic line depth dif- 
ferences. We see that the match of the Mg II A4481 pro- 
file is satisfactory, but the wings of the He I A4471 line 
appear too high, probably due to low level emission in 
this transition. Adding a 20% disk flux contribution and 
renormalizing the spectrum yields a diluted version ( dot- 
ted line), which appears to be much too weak compared 
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Fig. 6. — Hertzsprung-RusscU diagram, includingthe solutions 
from our small rapid rotator grid for /3 CMi (Tab. [Jll. The num- 
bers on the left side give the stellar mass of the corresponding 
evolution track, while the numbers at the bottom gives the evolu - 
tionary age correspond ing to the shown isoc hron es IYi et al.l l|2003l 'l: 
IDemarque et all II2004I '). As outlined in Sect. 14. ll only models A, B, 
and D satisfy the observational constraints on the stellar effective 
temperature. Furthermore, model D satisfies best our inclination 
obtained with interferometry (38.5 ± 1°). 



with the observed spectrum. This comparison suggests 
that the disk continuum emission in the blue range is 
negligible, and thus we will ignore any disk contribution 
to the optical flux in the following discussion. 

With the derived polar radius and bolometric lumi- 
nosity, we are able to locate the stellar position on the 
HR diagram, af ter correcting for the rotational effect 
(|Che et al.ll20Tll) . The stellar mass from the HR diagram 
is used as initial value for the next iteration step, until 
convergence between the assumed mass and the mass es- 
timated from the HR diagram is reached. 

For each parameter combination, our model provides 
the rotation period, apparent luminosity L'^pp, and ap- 
parent effective temperature T^^^, as listed in Tab. [2J 
Comparing the model effective temperature and luminos- 
ity with the observational con straints for /3 ClV Ii (Toff = 
12,050 K; L = 195 ± 60 L^: ISaio et alll2007l) . we find 
good agreement for models A, B, and D. Due to its con- 
sistency with the inclination and effective temperature 
(Sects. and H75)) . we favor model D, suggesting i ~ 
40°, w/wcrit = 0.99, Tpoio = 14, 550, w sin i = 244 kms^^, 
T^^^ = 12 871 K, and a polar and equatorial radius of 
i?poic = 2.95 Rq = 0.26 mas and i?cq = 4.10 Rq = 
0.36 mas, respectively. We plot the values corresponding 
to our model grid in the HR diagram shown in Fig. E] 
and co mpar e it to the evolut i onary tracks from lYi et alT 
(|2003f ) and IDemarque et al.l (|2004[ ). yielding an evolu- 
tionary age of 0.1 — 0.15 Gyr. 

Obviously, an important input parameter for our 
modeling procedure is the rotation velocity wsini, for 
which we use an av erage value from the literature. 
iTownsend et alj (|2004D argued that most v sin i measure- 
ments on Be stars might systematically underestimate 
the true projected rotation value due to the effect of 
gravity darkening. In order to test this scenario, we have 
artificially increased the measured u sin i- value by 10% 
and repeated our modeling procedure. Wc find that the 
increase in v sin i results in a significant increase in the 
stellar mass, apparent effective temperature, and appar- 
ent luminosity, making the model prediction much less 
consistent with the observed values. Therefore, wc sug- 
gest that in the case of /3 ClVIi, the v sin i-mcasurements 
are not significantly biased by gravity darkening, which 
is likely a result of the intermediate inclination angle 
of /3CMi (the effect would be stronger for an equator- 
on viewing angle) and the us e of our modera te gravity 
darkening law (Teff cx IChe et al.l [ 20ll . which is 

likely more relatistic than the classical v on Zeipel darken- 
ing coefficient (Tofi oc g^ff ^) adopted by ITownsend et all 
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(|2004f ). Clearly, the accurate determination of vaini re- 
mains a fundamental problem for constraining the stellar 
parameters of near-critical rotating stars. Therefore, we 
are currently working on incorporating the computation 
of model spectra in our modelling procedure, which will 
then be fitted to observed spectra together with the pho- 
tometric and intcrferometric constraints (Che et al., in 
prep.). 

The modeling provides spectra as well as model images 
in the continuum emission and in hydrogen photospheric 
absorption lines, which we will use as a representation 
of the stellar brightness distribution for our modeling of 
the disk in the following sections. 

4.2. Continuum disk geometry 

With projected baseline lengths up to B = 314 m. 
our CHARA/MIRC intcrferometric observations allow 
us to constrain the disk geometry in eight spectral 
channels in the _ff-band with an effective resolution 
of X/2B = 0.5 mas. A commonly applied model- 
ing approach is to approximate the disk emission with 
a Gaussian intensity profile, which is superposed on 
a uniform disk representing the stellar surface (e.g. 



Tvcner et al1[2005t iMeilland et al.l[200llPott et ai1[2010l : 



Schaefer et al.l 120101) . A caveat of this approach is that 
a significant part of the disk emission is concentrated 
at radii r < R^,, and therefore wrongly attributed to the 
disk instead of the stellar flux. Accordingly, these models 
will systematically overestimate the disk emission with 
respect to the stellar emission (i^disk/^*)j resulting in in- 
consistencies with SED fitting results. In order to avoid 
these biases, we employ an elliptical Gaussian model, 
where the radial intensity profile at r > is given by 
a half-Gaussian, resulting in a proper estimation of the 
^disk/-F* ratio. As alternative model, we considered an 
inclined ring model, where the radial intensity profile is 
given by a Gaussian centered at radius a from the star, 
with a fixed fractional width of 25% (see iMonnier et al.l 
l2006al for details). Based on radiative transfer simul- 
tations of the brightnes s distribution in classical Be star 
disks (e.g. 'Waters' 1986), we consider the Gaussian model 
a better representation of the expected emission profile in 
a Be disk, while the ring model is more sui ted for compar- 
ison w ith results from the literature (e.g. IMeilland et al.l 
[20M) . 

For all models, the stellar photosphere is represented 
by the aforementioned rapid rotator model (Sect. 14. ip 
and the position angle 6 of the stellar equator and the 
disk major axis are aligned. Besides the disk major axis 
a (ring major radius or Gaussian half-width half max- 
imum), the PA 9, and the inclination i, we also treat 
the flux ratio between the stellar and disk component 
(-Fdisk/-F*)H as a free parameter. The model was fitted to 
the MIRC squared visibilities, closure phases, and triple 
amplitudes using a least-square fitting procedure, result- 
ing in the best-fit model shown in Fig. [T] The best-fit 
model parameters are displayed in Tab. [3l including la- 
errors, which we have determined using a boot strapping 
procedure. We also tested whether the agreement can 
be imp roved using a skewed ring model (jMonnier et al.l 
l2006al) . but find that introducing any disk asymmetry 
only marginally improves the fit. Our model includes 
minor asymmetries due to the brightened polar region 
in our rapid rotator model, which results in small clo- 
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Fig. 7. — Left: Continuum visibilities measured on /3CMi with CHARA/MIRC in the H-hund (red points) and CHARA/CLIMB in the 
A'-band (blue points). In order to show the position-angle dependence of the visibility function, we have binned the data by the projected 
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dotted line) arc shown, with bot h mo dels yielding nearly indistinguishable visibility profiles. For ii'-band, wc show the best-fit Gaussian 
model (blue dashed line; see Scct. l4T2ll . Right: Intensity distribution corresponding to our best-fit ring (top) and Gaussian model (bottom). 



sure phases (< 1-2°), consistent with the measurement 
(Fig.©. 

The determined disk position angles of 139.2° (ring 
model) and 139.3° (Gaussian model) are in excellent 
agreement with the gas disk rotation axis determined 
with our VLTI/AMBER photocenter analysis (140.0 ± 
1.7°). 

Using the detailed information about the i7-band ge- 
ometry obtained with our extensive MIRC data set, we 
then fitted our Gaussian model to the CLIMB data 
to determine the ivT-band continuum geometry. Given 
the lower amount of observational constraints for the 
/C-band, we treated only the disk-to-star flux ratio 
{Fdisk/ Fh)k as a free parameters and kept the remaining 
parameters fixed. The resulting best-fit visibility curves 
are shown in Fig. [71 With (Fdisk/^;)^ = 0.25t^;J^^, we 
do not find any significant deviations between the K- 
band and iJ-band flux ratio. 

4.3. Investigating the differential phase signatures of 
the stellar rotation in the Brj absorption line 

With the currently achievable accuracy, DP measure- 
ments can already reveal photocenter displacements two 
to three orders smaller than the formal angular resolution 
iX/2B). For instance, our 2009-12-31 AMBER observa- 
tions exhibit a DP accuracy of ~ 0.8° (standard devia- 
tion over all continuum spectral channels), corresponding 
to a photocenter displacement of 8 micro-arcsecond 
(0.008 mas) on the employed 115 m baseline. Given that 
this is much smaller than the equatorial stellar diameter 
(~ 660 micro-arcsecond), it is important to investigate 
whether the measured DP signatures might also contain 
contributions from the photospheric Br7 absorption line. 



which is tracing the stellar rotation and is underlying the 
Br7-line emission. 

In order to simulate the effect of stellar rotation on our 
interferometric observables, we employ our rapid rotator 
code and compute the stellar surface brightness distribu- 
tion around the Br7-line with a similar resolution as our 
AMBER observations (Fig. [51 top). For this, we align 
the stellar rotation axis with the measured rotation axis 
(9 = 227.5°, Sect. [3|). Br7 absorption on one side of 
the photosphere will shift the photocenter towards the 
opposite direction. Taking this into account, we adjust 
the orientation of our model photosphere, so that the 
Br7 absorption in the blue-shifted line wing matches the 
measured direction of the photocenter displacement at 
red-shifted wavelengths. In addition to the photospheric 
emission, we include disk emission in our model, assum- 
ing {Fdisk/F*)K = 0.25, as determined in Sect. 14.21 

From the model, we compute the corresponding DPs 
and find signatures (j> < 0.5° on all baselines (Fig. [51 bot- 
tom). This result reflects the fact that the fraction of the 
photospheric absorption to the total K-hand emission is 
rather small (< 8%) and that the absorption appears 
as a rather diffuse structure on the extended stellar sur- 
face. We conclude that the photocenter signatures due 
to stellar rotation (in photospheric absorption lines) are 
a minor effect compared to the rotation signatures of the 
disk (in line emission) , but nevertheless result in DP sig- 
natures comparable to the currently achievable DP accu- 
racy and are therefore important, in particular for more 
detailed future studies. 

4.4. Constraints on the disk kinematics from the 
position- velocity diagram 



10 




0.2 0.0 -0.2 0.2 0.0 -0.2 0.2 0.0 -0.2 0.2 0.0 -0.2 0.2 0.0 -0.2 

dRA [mas] dRA [mas] dRA [mds] dRA [mds] dRA [mds] 



T 1 1 1 1 1 1 1 1 1 1 n 1 1 r 




I I I J I I I I I I I J I I I 

2.162 2.163 2.164 2.165 2.166 2.167 2.168 2.169 2.162 2.163 2.164 2.165 2.166 2.167 2.168 2.169 

Wavelength [p,m] Wavelength [)a,m] 

Fig. 8. — Top: Model images from our rapid rotator model for /3CMi for some representative wavelengths (Sect. [J^Sj. Bottom: Model 
differential phases (including the stellar photosphere only) computed for the baselines employed by our $ CMi VLTI measurements from 
2009-12-31 (left) and 2010-04-23 (right). 



Position-velocity (p — v) diagrams provide a powerful 
tool for the interpretation of a disk velocity field and are 
commonly employed in radio astronomy to derive the 
rotation profile of circumnuclear galactic disks or proto- 
stellar disks using, for instan ce, maser (e.g. lMivoshi et aTl 
ll995 HPestalozzi et a .l."2009') or molecular emission trac- 
ers (e.g. iSofue fcRl ibin 2001; Isella et al. 2003). Using 
our continuum-corrected photocenter displacement mea- 
surements, we can construct an equivalent diagram from 
our VLTI interferometric data by measuring the length of 
the continuum-corrected photocenter displacement vec- 
tors projected on the disk plane 9 = 140.0°. Performing 
this projection on the disk plane also allows us to avoid 
a potential bias due to opacity effects, since these effects 
would move the photocenter only perpendicular to the 
disk plane (Sect. [3]). The resulting p — w diagram (Fig. [3]) 
shows a symmetric rotation curve, which we interpret 
using the simpl e model of a thin Keplerian-rotating disk 
(|Weintroub et a l. 2008). The disk extends fr om an in- 
ner (i?in) to an outer (i?out) radius and the line-of-sight 
(LOS) velocity v of the emission element located at ra- 
dius r and at angle in the disk plane is given by 



Vkcp{r,d) = ^ —-^ sinii} ■ cosi, (4) 

where is measured against the LOS and G is the 
gravitational constant. Considering only the emission 
from the outer-most disk annulus (i.e. at i?out and "d = 
—7t...+tt) will result in a straight line in the p—v diagram 
(line A-B in Fig. [9]). Adding the emission from the re- 
maining disk annuli (i?in to i?out) will result in a charac- 
teristic "bowtie" -shaped filled region in the p-v diagram, 
such as commonly observed in CO ima ging observations 
(|Sofue fc Rubinll2001l:l!seria et al.l[200l . However, astro- 
metric observations (such as radio masers or photocenter 



displacements) trace the light barycenter of each annu- 
lus, corresponding t o the B-C and A-D cur ves in Fig. [9] 
(for more details see IWeintroub eTaI|[2008l) . We fit this 
simple model to the (3 CMi p — v diagram assuming a 
fixed inclination of i = 40° (as determined in Sect. 14. 2 1) . 
which yields best agreement for M^, = 3.6 ± 0.3 Mq, 
Rout = 2.4 ± 0.2 mas, and Rm < 1.0 mas (sohd curve in 
Fig.©. 

4.5. Detailed modeling of the disk kinematics 

By reducing the information content to purely astro- 
metric data, the p—v diagram analysis method presented 
in the last section provides a very intuitive and power- 
ful method to constrain the gas velocity field in the ob- 
served line tracer. In this section, we will use a more 
sophisticated model in order to make full use of the rich 
information contained in our spectro-interferometric ob- 
servations, including spectra, wavelength-differential vis- 
ibilities, DPs, and CPs. The Br7-emitting gas is assumed 
to be located in a thin disk plane, which is justified due 
to the expected small openin g angle of Be star disks (e.g. 
iBiorkman fc Cassinelli|[r993[ ) and the intermediate incli- 
nation angle under which /3CMi is observed (Sect. [4?2|) . 

The gas kinematics is parameterized with a rotation 
profile \v{r)\ = Acp(i?rof) • {r / R^cfY , where /3 = -1 for 
rotation with constant angular momentum, /? — —0.5 
for Keplerian rotation, /3 = for consta nt rotation , 
and 13 = +1 for solid body rotation (e.g. iSted 119961) . 
Acp = |i^(i?ref)|/|?4ep(^ref)| IS the Orbital vclocity at the 
reference radius i?i.ef = 1 AU, expressed as fraction of 
the Keplerian velocity \vkep{r)\ = (GM^/r)~^/^. In the 
model, the line emission extends from the stellar surface 
to an outer radius i?out with a radial power-law inten- 
sity profile Ii{r) oc r'. To include thermal line broad- 
ening in our kinematic model, we adopt a constant ra- 
dial gas temperature of 0.6 • T^s (=7860 K for /3CMi), 



Gas distribution, kinematics, and excitation structure in the disks around /3 CMi and C. Tau 



11 




200 



100 



> 

CO 

Q 




-100 



-200 



-1 1 
Position [mas] 

Fig. 9. — Left: Position-velocity diagram, derived from the continuum-corrected photocenter displacement vectors in Fig.\3\{bottom, left). 
The color of the data points matches the color-coding of the wavelength channe ls in Fig. [3] The model curves show the Kcplerian rotation 
profile for three stellar masses with = 0.68 mas, and Rout = 2.4 mas (Sect. BTill . Right: In o rder to i llustrate the interpretation of the 
p — ii-diagram, we show the Br7-line model images from our best-fit kinematical model f Sect. [431 Fig. llO|l for some representative velocities 
and mark the corresponding photocenter displacement (corresponding to the position offset plotted on the abscissa in the diagram) in 
each channel map with a blue cross. The yellow vertical line marks the position of the star, which corresponds to the photocenter of the 
continuum emission and the zero position in the p — v diagram. 



as suggested by the radiat ive transfer niodehng from 
iCarciofi fc BiorkmanI (|2006[) . 

We include both hne and continuum emission in our 
model and compute the interferometric observables for 
our given VLTI array configurations and the covered 
wavelength channels. The wavelength-dependent model 
visibilities and phases are then fitted to the VLTI inter- 
ferometric data using a redu ced Xr goodness-of-fit esti- 
mator (see iKraus et al.l 120091 for a definition), including 
our flux, visibility, DP, and CP constraints. 

For the continuum emission, we assume the geometry 
determined with our CHARA observations (Sect. 14. 2|) 
and use this model to renormalize the AMBER contin- 
uum visibilities. In order to incorporate the underly- 
ing photospheric Br7 absorption, we include the photo- 
sphere model discussed in Sect. 14.31 although, as dis- 
cussed above, the influence on the differential phase 
(< 0.5° on all baselines) is still within the measurement 
uncertainities. 

The disk position angle is fixed to 6* = 140.0°, as de- 
termined by our model-independent photocenter analysis 
(Sect. [3]). The remaining six parameters in our modeling 
are the outer disk radius -Rout, the inclination i, the stel- 
lar mass M^,, the radial intensity power-law index q, the 
disk rotation index /3, and the velocity at the reference 
radius i?rcf (expressed as fraction of the Keplerian veloc- 
ity, /kep)- In a first step, we fix the velocity profile to 
Keplerian rotation (/3 = —0.5, /kop = 1) and vary the re- 
maining four parameters systematically on a grid, yield- 
ing the best-fit model shown in Fig. [101 The correspond- 
ing parameters and uncertainties are listed in column 4 
of Tab. 131 In a second step, we test also non-Keplerian 
velocity fields, yielding the best-fit values listed in col- 
umn 5 (Fig. [T2l) . The resulting x^-surfaces are shown in 



Fig.[TTland the uncertainties have been derived using the 
bootstrapping technique. 

Our model fits show that the intriging phase inversion 
observed at the line center on our longest VLTI base- 
lines (Fig. [ini 3rd row) can be explained with the phase 
jumps in the Fourier phase crossing a visibility null. 
These phase jumps appear in the same spectral chan- 
nels where we measure visibility minima in the W^-shaped 
visibility profile, indicating that the visibility function of 
the pure line-emitting geometry transits here from the 
first to the second visibility lobe (Fig. [101 ^'^'^ and 3rd 
row). At the same time, the continuum emission is only 
marginally resolved, which results in the measured com- 
posite line+continuum visibility with a rather high con- 
trast of ^ 0.7. The phase jumps and visibility mini- 
mums are a basic property of the Fourier transform of 
strongly resolved objects and are reproduced very natu- 
rally and without finetuning from our kinematical mod- 
cling. The refore, our result s do not support the idea 
outlined by iStefl et all (|2011[ ) that these features might 
indicate secondary dynamical effects or polar mass out- 
flows. 

As bcst-flt value for the radial intensity index, we yield 
q = — 1.6±0.2. Isothe rmal viscous decretion disk rn odels, 
such as discussed in iBiorkman fc Carciofil ()2005() , pre- 
dict a radial disk surface density law E(r) oc r~^, cor- 
responding to q ~ —2 in the optically thin case. More 
realistic non-LTE disk models suggest a more shallow 
surface density profile (S(r) cx r"^ - "^) due to a steep 
tem perature drop in the i nner f ew stellar radii of the 
disk ICarciofi fc Biorkma^ (|2008f ). Therefore, we con- 
clude that our measured intensity profile is well consis- 
tent with these models. 

The most significant deviation of our simple kincmat- 
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ical model from the data is in the precise shape of the 
Br7-spectrum. For instance, the measured spectra show 
a weak emission component in the line center (Fig. llOp . 
which is not reproduced by the model and which might be 
related cither to the presence of a uniformly distributed 
low- velocity gas halo or, more likely, to radiative transfer 
effects. 

Another deviation between the model and the observa- 
tion concerns the weak V/ i?-asymmetry which can be ob- 
served in our VLTI/ AMBER Br7 spectra at both epochs 
(Fig. 131 top row) and which is not reproduced by our ax- 
ialsymmetric kinematical model. Within the 113 days 
covered by our AMBER observations, no changes in 
the V/R asymmetry could be observed, which is con- 
sistent with the co nclusion by iTvcner et al.l ()2005D and 
IJones et al.l (|2011D that the Ha profile does not show 
significant long-term variability. It is interesting to com- 
pare the average photocenter displacement between the 
blue- and red-shifted line emission (2.5 mas) with the 
characteristic size o f the Ha-emitting region (2.13 mas, 
iTvcner et al.l 120051 ) . which suggests that Br7 emerges 
from a similar spatial region as Ha and a much more ex- 
tended region than the H- and A'-band continuum emis- 
sion (-RBr7 ^ Rna > -^cont. Fig. , in agreement with 
the prediction from lCarciofil ()2010[ r However, a detailed 
comparison is difficult due to the non-simultaneity of the 
different observations. 

5. DISCUSSION ON C TAURl 

Tau is a particularly well-studied classical Be star 
which shows a cyclic variability in the flux ratio of 
the violet- and red-shifted wing of the double-peaked 



Ha emission line (jRivinius et al.ll2006l: IStefl et al.ll2009[) . 

These V/R variations exhibit a period of ^1429 days 
and are generally attributed to a global one-armed den- 
sity oscillation in a Keplerian (or nearly Keplerian) disk 
(IOkazakilll9"9l . 

A kinematic model for the Br7-line of some ear- 
lier AMBER MR observations on CTau was presented 
by IStefl et all (l200i) and ICarciofi et all ([20091) and 
suggested a Keplerian rotation profile. Our spectro- 
interferometric observations on C Tau provide a Br7 mea- 
surement at a new epoch in the Ha V/ R-phasc and cover, 
for the first time, also the wavelength region around the 
hydrogen Pfund lines. Besides our Br7 and Pfund series 
data, we include Ha-sizes from the lit erature, nam e ly the 
GET photocenter measurement by IVakili et al.l ()1998l 
^ 7 R^ or - 1.3 3 mas) and the M a rk 111 and N POl 
measurement s bv [ Quirrenbach et al] (|1994 I1997D and 
ITvcner et al.l (j2004l ). It is important to note that the 
Mark HI and NPOI results were based on visibility am- 
plitudes instead of differential phases. Also, a direct com- 
parison with these earlier observations is complicated by 
the known V/R variability of CTau, although we note 
that our measurement at phase 0.484 is reasonably close 
to the Ha-measure ment at phase 0.577, presented by 
ITvcner et all ()20Q4D . We convert the Gaussian FWHM 
derived by these studies to photocenter displacements 
by computing the centroid of the corresponding Gaus- 
sian brightness distributions. We find that Br7 emerges 
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is in good agr eement with the disk precession law proposed by 
ISchaefer et al.l 
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from a similar spatial region than Ha. but a more ex- 
tended region than the near-infrared continuum emis- 
sion, as already found for /3CMi fScct. I4.5p . The Pfund 
emission originates from intermediate stellocentric radii 

(^Ha ^Br7 > -Rpf > -Rcont; Fig. g]). 

5.1. Signatures of the known one-armed oscillation 

Both in Br7 and in the Pfl4-Pf22 lines, we detect a 
double-peaked line profile and clear rotation signatures 
in the differential phases. All line transitions exhibit a 
photocenter displacement with a stronger amplitude in 
the south-eastern (red-shifted) than in the north-western 
(blue-shifted) lobe (Fig. 21 middle panel). Such an asym- 
metric displacement is consistent with the presence of 
a one-armed oscilla tion in the disk (jStefl et al.l 120091 : 
Carciofi et ahl [20091) CHARA/MIRC observations by 



Schaefer et al.l (|2010f ) showed that the one-armed oscilla- 



tion pattern can also be observed as an asymmetry in the 
H-hand continuum emission. Using multi-epoch data, 
they find that the position angle of the asymmetry is cor- 
related with the spectroscopic variability and precesses 
around the star with the Ha V/R period. Our observa- 
tion (2010-01-01) adds an additional epoch for this analy- 
sis at the V/R phase 0.484 (assuming max imum phase at 
JD= 2454505.0 and a period of 1429 davs. lSchaefer et al.1 
l2010f ). We find that the spiral density maximum is lo- 
cated towards the south-east of the central star and the 
measured PA is cons istent with the sinusoidal PA mod- 
ulation proposed bv ISchaefer et al.l (|2010l ). as shown in 
Fig. [131 

The reported polariza, tion angl e for CTau (32.9 — 3 3.1°, 
IQuirrenbachiraII[T997t 35 ± 4° . [Ghosh et al.l[l99l 32- 
33.5°. IMcDavidlll999ir is in excellent agreement with the 
rotation axis position angle (35.8 ± 2.0) found by our 
spcctro-interferometric observations. 

5.2. Probing the disk excitation structure using 
multi-transition spectro-interferometry 

In order to test whether the measured differences in 
the stellocentric emitting radius of the Br7 and Pfund 
transitions arc consistent with the expected excitation 
structure in the disk, we construct a simple radiative 
transfer model assuming local thermodynamic equilib- 
rium (LTE). We assume that the line-emitting mate- 
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Fig. 14. — Radial intensity profile, as comp uted with our LTE 
model for different line transitions (Sect. 15.21 1. The intensity has 
been weighted by the emitting area and normalized to the peak 
intensity. 



rial is located in an equatorial disk which extends out- 
wards from the stellar radius with a constant vertical 
density per unit volume and a half-opening angle of 
8 = 5°. Then, we integrate for each radius r the op- 
tical depth in vertical direction, assuming hydrogen 
under LTE conditions and an isothermal temperature 
distributions. The number density for the different ex- 
citation levels and ionization stages is computed using 
the Saha and Boltzmann equation assumi ng a Gaussian 
line p rofile with thermal line broadening (jWilson et al.l 
|2009( ). The Einstein coefficients for the different hy- 
drogen transition s are estimated using the seri es expan- 
sion published by |Omidw^ScAUiiteB ([llll). The ra- 
dial temperature and surface density profile are param- 
eteriz ed with T(r) = Tcff{r/R^)-°-^ (jStee fc de Arauiol 
[1991 and Sfr) = Snfr/i?.)-^ with Sq = 2.1 gcm'^ 



( Carciofi et al.|[2009l ). For the stellar temperature, equa- 
torial stellar radius, and distance, we assume T^s = 
19 370 K, i?^ = 7.7 i?0, and d = 126 pc, respectively 
ijCarciofi et al.l[2009l ). 

Using the radiative transfer equation for LTE condi- 
tions {Iu{t) = Bi,{T){l — e~'^^'^^^)), we compute the emit- 
ted intensity per unit area as function of radius r in 
the disk, where is the Planck spectrum for temper- 
ature T at frequency v. From the radial intensity pro- 
files (Fig. [14]) we compute the centroid of the brightness 
distribution, which is then compared to the stellocentric 
emission radius measured in the different line transitions. 
Given that our model assumes a simplified vertical den- 
sity structure and does not include inclination effects, 
we do not aim to match the absolute sizes of the emit- 
ting region in all line transitions, but focus instead on 
the relative sizes. For this step, we normalize both the 
measured and the model photocenter offsets relative to 
Br7. The comparison between the Ha, Br7, and Pfl4-17 
relative sizes and our LTE model is shown in Fig. [T51 
We find that we can reproduce important observational 
features, in particular that 
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Fig. 15. — Origin of the Br7 and Pfund line emission from Tau, 
as measured by the length of the photocenter displacement vector, 
averaged over all velocity channels for a given line transition (the 
error bars represent the standard deviation in the different veloc- 
ity channels) . In addition, we include the Ha-line photocenter 
displacement byjVakili et al. (1998, blac k dat a point) and the Hct 
size e st imates from Q uirrenbach et al] II1994I ). IQuirrenbach et al.l 
119971) . ITvcner et aiTp004| ), which we converted from Gaussian 
FWHM to the corresponding photocenters (grey data point). The 
dashe d lin e shows the prediction from our LTE excitation model 
fSect. [53)l . 
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a) Br7 originates from a similar spatial region in the 
disk as Ha {RBr-y Riia)- This result is also in 
agreement with the predictions from more sophis- 
ticated non-LTE r adiative t ransfe r computations, 
such as made by ICarciofil (|2010[ ). In order to 
better characterize the differences between the 
Br7 and Ha-emitting region, contemporaneous 
observations in these two wavelength bands will 
be required. 

b) Br7 originates from a more extended region than 
the Pfund lines (i?Br7 > ^Pf)- Computing the 
weighted average of the measurements in the indi- 
vidual transitions, we yield i?pf/i?Br7 = 0.48±0.12. 
For the Pfund lines, we arc not aware of earlier pre- 
dictions obtained with non-LTE radiative transfer 
codes. 

Our observati onal results confirm the finding from 
iPott et al] (|201Clf ) , which found -Rbi^ > Rpf on the classi- 
cal Be star 48 Lib, and consolidates their suggestion that 
the measured size differences can already be explained 
with the expected optical depth differences between these 
line transitions. We encourage theoreticans and model- 
ers to employ their more sophisticated non-LTE radiative 
transfer codes in order to test the influence of inclina- 
tion and non-LTE effects, althou gh, based on the results 
from llwamatsu fc Hiratal (|2008l e.g. Fig. 3), we expect 
no significant departure from LTE for the inner disk re- 
gions and the high transitions traced by our observations. 
Future multi-transition spectro-interferometric observa- 
tions with improved uu-coverage might also measure the 
precise radial intensity profile in a model-independent 
fashion. Together with sophisticated radiative transfer 
simulations, these observations will reveal the excitation 
structure of the disk and constrain parameters such as 
the temperature profile and the vertical disk structure, 
which arc currently difficult to access. 

6. CONCLUSIONS 

Using CHARA and VLTI interferometry, our study 
combined high angular resolution (with baseline lengths 
up to 330 m) with kinematical information obtained at 
high spectral dispersion, yielding direct observational 
constraints on the gas distribution, excitation structure, 
and kinematics of the disks around two classical Be stars. 
Using a model-independent photocenter analysis method 
we derived the disk rotation axis for the prototypical ob- 
jects /3CMi and ^Tau and spatially and spectrally re- 
solved the disk rotation profile on scales of a few stellar 
radii. For both objects, we find that the determined 
gas rotation plane agrees well with the orientation of 
the continuum-emitting disk as resolved by CHARA, al- 
though there is also clear evidence for substructure in 
the disk around ( Tau, revealing a one-armed oscillation, 
as indicated by different displacement amplitudes in the 
blue- and red-shifted line wing. 

Using our data set on /3 CMi we constructed a position- 
velocity diagram, which can be interpreted using the 
well-established procedures from radio interferometry, 
but probes the milli-arcsecond scale position displace- 
ments resulting from the rotating disk around this Be 
star. From our kinematical constraints, we derive the 
dynamical mass of the central star to 3.5 ± 0.2, which 



is in excellent agreem ent with earlier spectroscopic stud- 
ies (|Saio et al.l 120071 ). The inclination of the system is 
38.5 ± 1°, as determined with our CHARA continuum 
and VLTI line observations. As shown with our de- 
tailed kinematical modeling, the rotation law is Keple- 
rian (/? = —0.5 ± 0.1) and we do not have to include 
an expanding velocity component in order to explain our 
data. Furthermore, our kinematical model allowed us to 
identify the origin of the phase inversion, which has now 
been observed in the differential phases in five out of eight 
Be stars. These phase jumps correspond to the transition 
from the first to the second visibility lobe, removing the 
necessity f or speculations be yond the canonical star-|-disk 
paradigm (jStefl et al.ll201lD . 

For ^ Tau, we obtained spectro-interferometric obser- 
vations covering simultaneously the Br7 and at least nine 
transitions from the Pfund line series. For all transi- 
tions, we detect a significantly stronger photocenter dis- 
placement in the red-shifted line wing than in the blue- 
shifted line wing, tracing the one-armed oscillation which 
has been deduced for the Tau disk before. Comparing 
the photocenter displacement in the different line tran- 
sitions, we find that the Pfund, Brackctt, and Balmer 
lines originate from different stellocentric emitting re- 
gions (i?cont < -Rpf < -RBr7 ^ Rlia) , which wc Can re- 
produce qualitatively with a simple LTE line radiative 
transfer model. More work, including non-LTE radia- 
tive transfer modeling, will be required in order to derive 
quantitative constraints. 

By detecting a purely Keplerian velocity field, our ob- 
servations are inconsistent with disk-formation mecha- 
nisms incorporating a strong outflowing velocity com- 
ponent, such as the wind compression scenario, which 
predicts a strong radial velocity co mponent compara- 
ble to the escape velocity (jBiorkman fc Cassinclli 1993). 
On the other hand, our kinematical constraints, as well 
as the measured hydrogen line intensity profiles (with a 
radial power law index q = — 1.6±0.2) are consistent 
with the pr edictions f rom Keplerian visco us de c retion 
disk models (|Lee et al.l [T991.). As sh own bvlKatol (I1983D 
and discussed in various reviews (e.g. lCarciofill201C ). vis- 
cous Keplerian disks are also able to produce one-armed 
density oscillations, such as detected for ^Tau. 

Considering that our kinematic constraints have been 
obtained using a very limited number of individual mea- 
surements (2 pointings on /3CMi, 1 pointing on ^Tau), 
our study also illustrates the high effectiveness achievable 
with spectro-interferometry, in particular if a very high 
spectral resolution is employed or several line transitions 
are observed. 
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